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Based on density functional theory (DFT) calculations including an on-site Hubbard U term we 
investigate the effect of substrate-induced strain on the properties of ferrimagnetic Fe203-FeTi03 
solid solutions and heterostructures. While the charge compensation mechanism through formation 
of a mixed Fe 2+ , Fe 3+ -contact layer is unaffected, strain can be used to tune the electronic prop- 
erties of the system, e.g. by changing the position of impurity levels in the band gap. Straining 
hematite/ilmenite films at the lateral parameters of Al2C>3(0001), commonly used as a substrate, 
is found to be energetically unfavorable as compared to films on Fe2O3(0001) or FeTiO3(0001)- 
substrates. 

PACS numbers: 73.20.-r,73.20.Hb,75.70.Cn,71.28.+d 



I. INTRODUCTION 

In the fabrication of ferromagnetic semiconductors 
for spintronics applications a lot of research focuses 
on the homogeneous doping of traditional or oxide 
semiconductors with 3d ionsi' 2 -' 3 -'^. However, the 
coupling between magnetic impurities and charge 
carriers is often too weak, leading to Curie temper- 
atures (Tc) way below room temperature (RT). On 
the other hand, materials like Fe2- :c Ti :r 03 exhibit 
intrinsic semiconducting and ferrimagnetic proper- 
ties, although the end members a-Fe2 03 and FeTi03 
are antiferromagnetic insulators with Tn = 948 and 
56 K, respectively. Besides applications in spintron- 
ics, this material is also discussed in paleomagnetism 
as a possible cause of anomalies in the Earth's mag- 
netic field, as well as for electronics devices (e.g. 
varistors) because it is a wide band gap semiconduc- 
tor that can be either n- or p-type depending on the 
doping concentration^. A Curie temperature above 
RT and a reduction of resistivity was observed in 
synthetic solid solutions with Ti concentrations up 
to 70%S£ Moreover, T c was found to increase upon 
annealing both in these samples and in thin epitax- 
ial films 8 . This behavior can be attributed to cation 
ordering phenomena related to a miscibility gap in 
the rather complex phase diagram of the system^. 

The origin of ferrimagnetic behavior remained 
unclear until recently. Both materials have a 
corundum(-related) structure (see Fig. [1]) with a 
stacking of 2Fe 3+ /30 2 ~ in hematite (space group 
R3c) and 2Fe 2+ /30 2 -/2Ti 4 +/30 2 ~ in ilmenite (i?3) 
along the [0001]-direction. Thus at an interface or in 
a solid solution (SS) charge is not compensated, if 
all ions preserved their bulk valence states. DFT 
calculations considering correlation effects within 
LDA+ttiS showed that the charge mismatch is ac- 
commodated by a mixed Fe 3+ , Fe 2+ contact layer at 



the interface^, providing first theoretical evidence 
for the lamellar magnetism hypothesis^. The Fe 2+ - 
ions at the interface give rise to uncompensated mo- 
ments and also to impurity states in the band gap. 

The incorporation of Ti in hematite^ 3 , (a = 5.04 
A, c = 13.75 A) introduces a substantial strain: the 
volume of the end member ilmenite^ (a = 5.18 A, 
c = 14.27 A) is 9.7% larger than the one of hematite. 
Indeed, lens-shaped dark contrasts around nanoscale 
hematite lamellae in an ilmenite host, imaged by 
transmission electron microscopy, indicate signifi- 
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FIG. 1: (Color online) Crystal structure of the 60-atom 
unit cell of Fe2-^Ti a: 03 for x — 0.33 with a layered (a) 
and more homogeneous arrangement of the Ti-cations 
with Ti in the same (b) and different (c) spin-sublattices. 
Oxygen, Fe and Ti are shown with light grey, red and 
black spheres respectively. Pink circles mark the Fe 2+ - 
positions, while the rest of the iron are Fe 3+ . The local 
magnetic moments at the cation sites and the total mag- 
netization of the system in fj,B are given in the right side 
and bottom of each configuration, respectively. 



cant strain fields^ 2 -. 

Epitaxial Fea-^Ti^Oa films^^^^^ 



are typ- 



ically grown on an Al2O3(0001)-substrate (a = 
4.76 A, c = 12.99 A) which introduces a substantial 
compressive strain of 5.8% and 8.8% compared to 
Fe 2 03 and FeTi03 and only rarely, a Cr 2 03-buffcr 
layer is used^ to reduce the lattice mismatch. 

Epitaxial strain can have a strong impact on the 
film properties, e.g. by tuning the magnetic inter- 
actions in magnetoelastic composites 20 , enhancing 
ferroelectricity 21,22 or even inducing orbital recon- 
structions^. The goal of the present study is to 
explore the effect of strain on the properties of 
Fe2- a; Ti a ;03. In particular we address its influ- 
ence on (i) the energetic stability and compensation 
mechanism as well as on (ii) the electronic, mag- 
netic and structural properties of the system. DFT 
calculations are performed on SS and layered config- 
urations with x = 0.17,0.33,0.50 and 0.66, strained 
laterally at the lattice parameters of AI2O3, Fe203, 
and FeTi0 3 . 



II. CALCULATIONAL DETAILS 

We use the all-electron full-potential linear aug- 
mented plane wave (FP-LAPW) method as im- 
plemented in the WIEN2K code^ and the gener- 
alized gradient approximation (GGA) 25 . Within 
LDA+U 10 U = 8.0 eV and J = 1.0 eV is applied 
to the Fe and Ti 3d states. These values were found 
to reproduce correctly the ground state of FeTiO^. 
The systems are simulated in a hexagonal unit cell 
with 60 atoms (Fig.QJ. Besides the layered configu- 
rations (cf. Fig.[T^,) more homogeneous distributions 
are generated by substituting 50% of Fe in a bilayer 
by Ti, as shown e.g. in Fig. [TJa-c. For further details 
on the calculation see (Ref^). 



III. RESULTS AND DISCUSSION 

The optimized c/a-ratio and volume (Fig. [2^i-b) 
show a linear increase with xn m in accordance with 
Vegard's law, similar to what was observed exper- 
imentally in synthetic hematite-ilmenite solid solu- 
tions 6 . Furthermore, for a given concentration both 
c/a and V are largely independent of the distribu- 
tion of Ti- impurities. The c/a-ratio of bulk FeTiC>3 
(2.76) is slightly larger than the one for a-Fe203 and 
AI2O3 (2.73). Due to the small tensile/compressive 
strain when using aFeTi0 3 /aFe20 3 the c/a-ratio of 
Fe 2 - x Ti x 3 is slightly reduced (-1.1 to -2.8 "^/in- 
creased (3.1-5.2 %), respectively. In contrast, due to 
the high compressive strain on an A^Os-substrate, 
c/a increases strongly by 14.7-16.6 % which corre- 
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FIG. 2: (Color online) (a) c/a-ratio (b) volume and (c) 
formation energy (eV/f.u) versus ilmenite concentration 
xn m for Fe2-:zTb;03 strained at the AI2O3 (red/dark 
grey), Fe2C>3 (grey) and FeTiOa (black) lateral lattice 
constants. Circles (triangles) denote compensation in- 
volving Ti 4+ (Ti 3+ ). Open/filled symbols refer to solid 
solutions (SS)/ layered configurations (L). Horizontal 
lines mark the bulk c/a ratio and volume of the end 
members and AI2O3. Red (dark grey) squares indicate 
experimental data from Takada et al— . 



sponds to c re i — 14.89 — 15.15 A. Nevertheless, the 
volume does not completely relax: The volume of 
the system strained at the Al2 03-lateral lattice con- 
stant is 6.8 % (10.2 %) smaller than when strained 
at a Fc2 o 3 (aFcTio 3 )- Tne volumes of Fe 2 - x ^x0 3 
strained at aF e Ti0 3 and aFe 2 3 lie between the ones 
for the end members Fe203 and FeTiOs. 

X-ray diffraction data for Fe2- a; TL r 03 films on 
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Al2O3(0001)2£>2£ indicate significant lateral strain 
relaxation: already in a 10 nm thick film a relaxes to 
the bulk value of FeTiC>3 with only a small change 
in c/a (see Fig. The c/a values and volumes ob- 
tained by Takada et aii^I are in good agreement with 
the DFT values of the systems strained at aF e Ti0 3 - 

Next we turn to the influence of strain on the en- 
ergetic stability. The formation energy with respect 
to the end members as a function of ccxi is shown 
in Fig. for the three different substrate lattice 
constants. For each Ti-concentration we have con- 
sidered several different cation arrangements, e.g. 
for x = 0.33 these include an ordered arrangement 
with an Fe layer sandwiched between two Ti layers 
(Fig. QJi) or solid solutions with Ti ions either in the 
same (Fig.njj) or different spin-sublattices (Fig.QT;). 
We find that compensation through Ti 4+ and dispro- 
portionation in Fe 2+ , Fe 3+ is more favorable over 
mechanisms involving Ti 3+ . Furthermore, the for- 
mation energy increases linearly with xn m . These 
features are independent of the substrate lattice pa- 
rameters. Systems strained laterally at ap C Tio 3 are 
more stable than the ones on aFc 2 o 3 ■ ln contrast, the 
formation energy of films strained at aAi 2 o 3 increases 
by 0.7 eV as compared to films on aFe 2 o 3 - This im- 
plies that the strong compressive strain is energet- 
ically unfavorable and gives a possible explanation 
why a lateral strain relaxation occurs in Fe2- a; Ti a ;03 
film o 26 i 27 . While for systems strained on hematite 
and ilmenite substrates layered arrangements (full 
symbols) are more favorable than homogeneous dis- 
tributions (open symbols), the trend is reversed for 
x = 0.33 and x — 0.66 on an Al 2 O3(0001)-substrate. 

Concerning the electronic properties of the hemo- 
ilmenite system, we have plotted in Fig. [3] the den- 
sity of states of a Ti-double layer in a hematite host 
(Fig. QJi), but similar behavior is observed for all 
studied systems. Upon Ti 4+ substitution, an iron- 
ion from the neighboring layer turns Fe 2+ , as ob- 
served also for isolated impurities by Velev et al~. 
The Fe 2+ 06 and the TiOe-octahedron are corner- 
(and not face-) sharing. The so formed Fe 2+ -ions 
in the contact layer have an impurity state of a\ g 
symmetry (d z 2) that is pinned at the Fermi level 
for systems strained at aF C2 o 3 and aFeTiCV Such 
a mid-gap state was recently reported from x-ray 
valence band photoemissior* 2 ^ and optical measure- 
ments^, although it was related to the low oxygen 
pressure during deposition. The main feature re- 
lated to strain is the change in band width: While 
for tensile strain at aF e Ti0 3 the bands are narrowed, 
for compressive strain at aAi 2 o 3 they are strongly 
broadened. This results in a reduction of the band 
gap (between the impurity state defining the Fermi 
level and the bottom of the conduction band) from 
1.90 eV for ap C Ti0 3 and 1-79 eV for aF C2 o 3 to 1.43 eV 
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FIG. 3: (Color online) Density of states of Fe1.67Tio.33O3 
containing two Ti-layers in a hematite host (shown in 
Fig.rjJ,): a) total; b) and c) projected of the 3d states of 
Fe 2+ at the interface and between the two Ti layers. The 
DOS of the system strained at the lateral lattice param- 
eters of Fe203, AI2O3, and FeTi03 is shown with a grey 
shaded area, red (dark grey), and black line, respectively. 



for aAi 2 o 3 - The corresponding values for x — 66% 
show the same trend but are smaller: 1.64 eV for 
aFeTi0 3 and 1.46 eV for aFe 2 o 3 to 0.78 eV for aAi 2 o 3 - 
The local magnetic moments and total magneti- 
zation for the three systems with x = 0.33 is dis- 
played in Fig. [1] Strain has only a small impact 
on the magnetic moments of Fe 2+ (~ 3.5 and 
Fe 3+ (~ 4.1 /xb) respectively which are reduced by 
less than 0.05 /xb at aAi 2 o 3 - The Fe 2+ -layer sand- 
wiched between two Ti-layers in Fig. [lji is only 
weakly coupled to the next Fe-layer (parallel and 
antiparallel orientation of the magnetic moments is 
nearly degenerate as in the ilmenite end member). 
Therefore, at temperatures above the Neel temper- 
ature of ilmenite, such layers will not contribute to 
the total magnetization. In contrast, Fe 2+ in the 
contact layer shows a strong antiferromagnetic cou- 
pling to the neighboring Fe-layer of the hematite 
host. These defect interface moments are respon- 
sible for the ferrimagnetic behavior of the system 
(M to t = 8.0 hb)- In solid solutions, Ti substitution 
in different spin-sublattices (e.g. in adjacent layers 
as shown in Fig.[T};), resulting in a zero net magne- 
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tization, is less favorable compared to substitution 
in the same spin-sublattice (Fig. [TJd), which maxi- 
mizes the total magnetization (M to t = —16.0 /xb). 
This trend promotes ferrimagnetic behavior in the 
system. 

IV. CONCLUSIONS 

Density functional theory calculations within 
GGA+t/ show that the charge compensation in 
hematite-ilmenite heterostructures and solid solu- 
tions takes place through a mixed Fe 2+ , Fe 3+ con- 
tact layer. This mechanism is robust with respect 
to substrate-induced strain. For Fe2C>3(0001) or 
FeTiO3(0001) substrates layered arrangements are 
more stable than solid solutions. However, the com- 
pressive strain at aAi 2 3 is likely to cause a stronger 
competition and even reverse the trend for x = 0.33 
and x — 0.66. The growth of epitaxial films on 



an A^Ovsubstrate is connected with a high energy 
cost. Therefore, in order to release strain such films 
may roughen or buckle in the first layers as recently 
reported by Popova et al£±. In contrast, the growth 
on lattice matched substrates or even substrates that 
produce a small tensile strain like FeTiO*3 is energet- 
ically favored. Our DFT results indicate that strain 
can have a strong impact on the structural and elec- 
tronic properties in the hematitc-ilmcnite system: 
e.g. by tuning the band width or the position of im- 
purity levels in the band gap and thus changing the 
concentration of spin-polarized carriers. 
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